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There has been a growing interest in realizing topologically nontrivial states of matter in band insulators, 
where a quantum Hall effect can appear as an intrinsic property of the band structure. While the on-going 
progress is under way with a number of directions, the possibility of realizing novel interaction-generated topo¬ 
logical phases, without the requirement of a nontrivial invariant encoded in single-particle wavefunction or band 
structure, can significantly extend the class of topological materials and is thus of great importance. Here, we 
show an interaction-driven topological phase emerging in an extended Bose-Hubbard model on kagome lattice, 
where the non-interacting band structure is topological trivial with zero Berry curvature in the Brillouin zone. 
By means of an unbiased state-of-the-art density-matrix renormalization group technique, we identify that the 
groundstate in a broad parameter region is equivalent to a bosonic fractional quantum Hall Laughlin state, based 
on the characterization of unverisal properties including groundstate degeneracy, edge excitations and anyonic 
quasiparticle statistics. Our work paves a way of finding interaction induced topological phase at the phase 
boundary of conventionally ordered solid phases. 


I. INTRODUCTION 

While fundamental particles in nature are either bosons 
or fermions, the emergent excitations in two-dimensional 
strongly-correlated systems may obey fractional or anyonic 
statistics'’^. The most famous example is the quasiparticles or 
quasiholes in fractional quantum Hall (FQH) effects in strong 
magnetic field, which are topological states of matter whose 
low-energy physics is governed by the Chem-Simon gauge 
theory. Interestingly, recent theoretical discoveries have re¬ 
vealed that FQH effects can also be realized in lattice systems 
in the absence of an external magnetic field^“^. Such a lat¬ 
tice realization of FQH phases is attributed to two key points: 
a nearly dispersionless single-particle energy band with non¬ 
trivial topology characterized by a nonzero Chern number and 
a strong many-body interaction comparing to the band width. 
The condition of nearly flat topological band is esentially im¬ 
portant, as the kinetic energy of particles can be quenched in 
such a topological band akin to the Landau level physics. The 
strong interaction also plays the vital role in stabilizing the 
FQH phases. In fact, without the interaction or interaction be¬ 
ing a subleading correction, the system is expected to be in a 
Fermi liquid like state at fractional fillings rather than forming 
a topological phase. 

Given above facts, a natural question that we address in 
this paper is whether a FQH phase is also possible in a lat¬ 
tice model with trivial non-interacting band. Actually, there 
have been a series of proposals along this direction^”^^, where 
the common wisdom is that the presence of strong interac¬ 
tions in a strongly frustrated system can give rise to a nonlocal 
complex bond order parameter and a spontaneous breaking 
of time-reversal symmetry (TRS) through flux attachment^. 
However, lacking of a theoretical method to predict the quan¬ 
tum phase for microscopic interacting systems, theoretical 
studies usually resort to different mean-fleld approximations, 
which often favor topological phases. As an example, Raghu 
et. showed that a quantum anomalous Hall effect can be 
dynamically generated in an extended Hubbard model on the 


honeycomb lattice. The similar idea has been applied to other 
lattice systems with a quadratic band crossing point, such as 
kagome^\ checkboard^^, diamond^^ and Lieb lattice^"^. How¬ 
ever, comprehensive numerical studies have been searching 
for true groundstates in different lattice systems and failed 
to And exotic topological phases predicted by the mean-fleld 
theories Here a crucial difficulty is that, instead of trig¬ 
gering the desired TRS spontaneously breaking, the strong 
interactions also tend to stabilize competing solid orders by 
breaking translantional symmetry. Thus, the simple concept 
of realizing interaction-induced FQH phases in topological 
trivial bands was illusive in realistic lattice models. 

Very recently, theoretical studies of extended kagome anti¬ 
ferromagnetic systems have discovered a particular class of 
spin liquids, the so-called chiral spin liquid^^’^^, with TRS 
spontaneously^^’^^ (or explicitly^^) broken, which shed lights 
on this elusive area: Long-ranged frustrated interactions may 
favor the FQH-like ground state near the boundary between 
ordered states^^. So far the existing examples are rare and 
ah occur at the half-flhing (half of spins are pointing up in 
z-direction) on kagome lattice, which may be attributed to 
quantum fluctuations near the non-coplanar spin ordered state 
(cuboc phase)^^”^^. Thus it is highly desired to search for the 
interaction-induced FQH phase beyond the half-flhing, which 
serves as the proof of the principle that TRS broken phase can 
emerge in more general conditions without a nearby cuboc 
phase. On the other hand, kagome-based magnetic systems 
have been widely studied under an external magnetic fleld^^, 
which can tune the spin systems into different magnetizations 
corresponding to hardcore boson systems at different flhings. 
The interesting candidate states have been established includ¬ 
ing the valence bond crystal state^^”^^ and the featureless Mott 
insulator^^ as possible groundstates. Beside these topological 
trivial phases, a Z 2 topological phase may survive in an easy- 
axis kagome system^^, which is currently under debate^^. It is 
theoretically proposed that FQH state can also emerge at 1/3 
fllling^^’^^, however, so far this possibility has not been estab¬ 
lished by controlled theoretical methods beyond mean-fleld 
approaches. Along this line, the existence of a topological 
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ordered phase at one-third filling remains open. 

In this paper, we study an Bose-Hubbard model on the 
kagome lattice in the hard-core limit: 


H = t ^ ^bl,br + H.C. - 1 - Vi ^ TlrTlr' 

(rr') (rr') 

+ 1^2 ^ rirflr' + V^3 ^ TlrTlr' (1) 


((rrO) 


(((rr'))) 


, where (br) creates (annihilates) a hard-core boson at site 
r. t = 1 is the nearest-neighbor hopping amplitude, , V 2 , V 3 
are the density-density repulsion strengths on first, second and 
third nearest neighbors, respectively. We focus on the bo¬ 
son filling number 1 /3 in this paper. This model can also be 
mapped onto the spin—1/2 XXZ model, allowing for an inter¬ 
pretation of our results in terms of both bosons and quantum 
spins. The energy band for hosting hardcore bosons is topo¬ 
logical trivial (with real hopping terms) with zero Berry cur¬ 
vature in the Brillouin zone. In order to study the ground state 
phase diagram in the {Vi, V 2 , V 3 } parameter space, we im¬ 
plement the density-matrix renomralization group (DMRG) 
algorithm on cylinder geometry combined with the exact di- 
agonalization (ED) on torus geometry (see Appendix A for 
computational details), both of which have been proven to be 
powerful and complementary tools for studying realistic mod¬ 
els containing arbitrary strong and frustrated interactions. 


II. RESULTS 
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FIG. 1 : Phase diagram of the extended Bose-Hubbard model on 
kagome lattice obtained by DMRG calculations on cylinder of cir¬ 
cumference Ly = 4 . (a)The phase diagram plotted in Vi and 
V2 = V3 parameter space. The shaded area is a coexistence re¬ 
gion. (b)The phase diagram plotted in V2 and V3 parameter space by 
setting 1/1=0. The contour plots of static density structure factor 
for: (c) charge density wave q = ( 0 , 0 ) phase, (d) strip phase and (e) 
FQH phase. The white dashed line shows the first Brillouin zone. 


near the classical phase boundary neighboring with different 
solid phases (Appendix B). 


A. Phase Diagram 

Our main findings are summarized in the phase diagram 
Fig. l(a-b). In the parameter region V 2 ~ V 3 and I/ 2 , V 3 > Vi 
(0 ^ Vi ^ 2.0), we find a robust FQH phase emerging with 
the TRS spontaneously breaking. The FQH phase is centered 
around the line V 2 ~ V 3 , which persists to IV 2 — V 3 I <0.1 
approximately as shown for Vi = 0 in Fig. 1(b). The FQH 
phase is characterized by a four-fold groundstate degeneracy 
on torus geometry, which arises from two sets of Laughlin 
z/ = 1/2 FQH states with opposite chiralities. In particular, 
the TRS spontaneous breaking can be inspected by identify¬ 
ing local circulating currents on the cylinder geometry, while 
the nature of bosonic Laughlin u = 1/2 state will be identi¬ 
fied by the edge excitation spectrum, fractional Chern number 
and the anyonic quasiparticle statistics as elaborated later. In¬ 
terestingly, we also show that the FQH liquid phase is neigh¬ 
boring with several solid phases which all respect TRS: a strip 
phase, a charge density wave q = ( 0 , 0 ) phase, and q = ( 0 , tt) 
phase with q as the ordering wave vectors. Compared to the 
distinctive Bragg peaks in the structure factor for solid phases 
(Fig. l(c-d)), the FQH phase shows a structureless feature 
(Fig. 1(e)). Finally, between the FQH phase and q = (0, 0) 
phase, there exists a narrow window for the coexistence of 
both FQH and charge density order (labeled by shaded area in 
Fig. 1(a)). By comparing our quantum phase diagram with the 
classical phase diagram, we find that the FQH phase is present 


B. Energy Spectrum and Doubled Topological Degeneracy 

Topological ordered states have characteristic groundstate 
degeneracy on compactified space (i.e. torus) while TRS 
spontaneous breaking topological phase has doubled topolog¬ 
ical degeneracy. To demonstrate this property in the interme¬ 
diate FQH region, we first investigate the low-energy spectra 
based on ED calculation. In Fig. 2, we show the scan of en¬ 
ergy spectra along the line Vi = 0.5 and varying V 2 (= Vs)- 
It is clear shown that there is a fourfold groundstate degen¬ 
eracy in the regime 0.8 < V 2 = V 3 < 1.5, which is sepa¬ 
rated by higher excited states by a robust spectrum gap. The 
fourfold degeneracy arises from two-fold topological degen¬ 
eracy for the V = 1/2 Laughlin state (full evidences will be 
shown below) and an additional factor of 2 from two-copies 
of states with the opposite chirality due to the TRS sponta¬ 
neously breaking. All four groundstates are located in mo¬ 
mentum sector k = ( 0 , 0 ), consistent with the expectation 
of the momentum folding rule for emerging FQH state with 
12 particles on 3 x 3 x 4 kagome lattice^’^. The low-energy 
spectrum gap for different system sizes is shown in Fig. 2(b), 
where we find that the energy gap between the fourth lowest 
energy state and the fifth one is robust against the increase of 
system sizes (Fig. 2(b)). This result indicates the emergent 
FQH phase may be robust at thermodynamic limit, which will 
be further confirmed by our larger system results based on 
DMRG. We also find that the energy gap is robust in the whole 
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FIG. 2: Energy spectra from ED calculation, (a) A scan across the 
reference line in Eig. 1 by setting Vi = 0.5 on the = 3 x 3 x 
4 = 36 sites cluster. The four-fold degeneracy of the EQH state 
is present around 0.8 < V2 = V3 < 1.5. Different colors and 
symbols correspond to differeent momentum sectors, (b) Energy gap 
for various system sizes Ns by setting Vi = 0.5, V 2 = V 3 = 1.2. 
(c)Contour plot of energy gap versus Vi and V2 = V3 on Ng = 36 
sites cluster. 


FQH regime while it drops to near zero at the phase boundary 
as illustrated in Fig. 2(c). 



FIG. 3: Local current Jij pattern of from DMRG calculation. 
The parameter is lA = 0.5, 1/2 = V 3 = 1.2 on = 4 cylinder 
(only show the five columns in the center). Width of bond is propor¬ 
tional to the current value (Jij is shown on the bond as a number) 
and arrows correspond to current directions. 


C. Time Reversal Symmetry Spontaneously Broken and Local 
Current 

Since the Hamiltonian (Eq. 1) preserves TRS, the emergent 
FQH state as the groundstate breaks TRS spontaneously. To 
investigate this mechanism, we move to larger systems with 
the cylinder geometry and obtain the groundstates 
by implementing DMRG calculation. By using different ran¬ 
dom initial wavefunction in DMRG, we obtain four different 
groundstates. Here we label topologically different ground- 
states by their chiral and anyonic nature on cylinder geome¬ 
try, where L (R) stands for “left-hand” (“right-hand”) chiral¬ 
ity, and 11 (s) stands for identity (semion) quasiparticle (see 
below).) A simple picture of the TRS broken can be obtained 
by measuring the local circulating currents in real-space. As 


shown in Fig. 3, the current Jij = 16^6^ | ) between 

two nearest-neighbor sites (i, j) forms loop structure and is in 
the clockwise direction in each triangle in the bulk of the sys¬ 
tem, which is refered as “left-hand” chirality. The emergent 
loop current is a direct demonstration of TRS spontaneously 
breaking for the state, which enables experimental detecting 
from local current measurements. This is a distinct properties 
different from the chiral spin liquid at half filling, where only 
the three spin chirality term is nonzero^^’^"^. We also find that 
the TRS partner hosts anti-clockwise loop current in 

each triangle. 


D. Fractionalization and Fractional Statistics for 
Quasiparticles 

To uncover the anyonic nature of groundstates in the in¬ 
termediate region, we investigate hallmark signatures of FQH 
state including characteristic excitation spectrum on the edge, 
fractional Chern number and quasiparticle braiding statistics 
in the bulk. All the evidences we obtain fully support the 
topological phase in the intermediate region is the emerging 
bosonic 0 = 1/2 Laughlin state. 

1. Entanglement Spectrum 

Firstly, we study the characteristic edge excitation with the 
help of entanglement spectrum (ES)^^, as partitioning a cylin¬ 
der into two halves manifests a “spatial” boundary. Fig. 4 
shows the ES for two of the groundstates with “left” 

chirality. The ES is grouped by the relative boson num¬ 
ber NNl of the half system and their relative momentum 
quantum number (relative to the total Ky of the high¬ 
est weight spectrum level) along the transverse direction (re¬ 
ferred to as y-direction). The leading ES of ^ displays 

the sequence of degeneracy pattern {1,1, 2,3, 5,...} in each 
NNl sector. Importantly, the edge mode counting rule agrees 
with the prediction of the free chiral boson in S'(7(2) 1 confor¬ 
mal field theory which describes the edge theory of Laughlin 
state^^. In addition. Fig. 4 also signals the chiral nature of the 
edge spectrum (ES increases as Ky varies from 0 to — 27r by 
a step dKy = —27r/I/^), which results from the TRS spon¬ 
taneously breaking. The other two groundstates with 

“right” chirality have opposite chirality {5Ky = 2iT/Ly) but 
the same degeneracy pattern in ES (not shown). 


2. Fractional Charge and Chern Number Quantization 

Secondly, we perform a numerical flux insertion simula¬ 
tion on cylinder systems^^’^^“^\ to determine the quantiza¬ 
tion of Hall transport and the topological Chern number of 
the ground state. This simulation realizes Laughlin gedanken 
experiment"^^"^"^, where a quantized charge will be pumped 
from one edge to the other edge by inserting a (7(1) charge 
flux in the hole of the cylinder for a quantum Hall state. As 
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FIG. 4: The low-lying ES of (left) and (right) on = 6 
cylinder. The ES is labeled by the relative boson number ANl = 
Nl — of left half cylinder in each tower (N^ is the boson num¬ 
ber of the state of reduced density matrix with the largest eigen¬ 
value). In each tower, the horizontal axis shows the relative mo¬ 
mentum AKy = Ky — Ky in the transverse direction of the corre¬ 
sponding eigenvectors {Ky is momentum of the state with the largest 
eigenvalue in each tower). The numbers below the red dots label the 
nearly degenerating pattern for the low-lying ES with different AKy. 
The black dashed line shows the entanglement gap in each momen¬ 
tum sector. 


shown in Fig. 5 (a), by threading a flux quantum 0 = 27r, 
adiabatically evolves into ). Further increasing flux 

up to Ttt will drive the system back to the |^^)- Interest¬ 
ingly, comparing the ES at = 0 and Ttt, the adiabatic flux 
insertion shifts the lowest level of ES from AN^ = 0 to 
ANl = —1, signaling a net charage transfer AQ = 1 (a 
unit charge) from the left edge to the right edge. In Eig. 5(b), 
the net charge transfer AQ is demonstrated, which is nearly 
quantized at AQ ^ 0.50 at 6> = 27r. Based on these observa¬ 
tions, we identify the bulk Chem number of the groundstate 
asC^ = Cs = 1/2, fully characterizing the obtained state as 
the Laughlin u = 1/2 state. 

Here we identify each groundstate hosts a fractional Chem 
number C = 1/2, which is similar to the Laughlin v = 1/2 
state in fractional Chern insulator"^^. The key difference is, the 
ground state in the fractional Chem insulator inherits the non¬ 
trivial topology from the non-interacting band. However, our 
extended Bose-Hubbard model (Eq. 1) has a Chem number 
C = 0 for the single-particle band. Thus, our model realizes 
an interaction-driven topological phase from a topologically 
trivial band structure"^^ resulting from the spontaneous TSR 
breaking. 


3. Modular Matrix 



FIG. 5: (a) The ES flow with inserting flux 0 in the hole of the cylin¬ 
der with Ly = 4. Starting from the identity ground state at 

^ = 0, the system evolves into the ground state by adiabati¬ 
cally threading a ^ = 27r flux. Eurther increasing flux up to An will 
drive the system back to the identity ground state (b) Charge 

transfer from one edge to the other edge on the cylinder geometry. 


ent quasiparticles. The Z//—matrix encodes the self-statistics 
of the quasiparticles, i.e. topological spin hi. Here we uti¬ 
lize the route of the “twist” overlap between the two topolog¬ 
ically degenerating groundstates and to constmct 
the modular S and U matrices^^”^^. 

The obtained results at Vi = 0.5, V 2 = Vs = 1.0 on = 4 
cylinder are 
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Indeed, the numerical obtained modular matrices are very 
close to the analytical prediction from SU (2)i Chem-Simons 


theory47-49: ^ 


1 1 
1 -1 


and = 


e 


1 0 
0 i 


From S— and U— matrices, we obtain full 


statistics of fractionalized quasiparticles: 1) There are two 
kinds of quasiparticles in total: Identity 11 and semion s, and 
the total quantum dimension is V = a/2; 2) The fusion mle of 
quasiparticles (that specifles how the quasiparticles combine 
and fuse)"^^: llX 5 = 5 , 5 Xs = ll;3)ll and 5 respectively 
has topological spin (the phase factor for the quasiparticle ob¬ 
tained during a self-rotation of 27r): = 0, hg = 1/4. 

These braiding statistics provides the strongest confirmation 
that the topological groundstate is equivalent to the bosonic 
FQH = 1/2 state. 


We further demonstrate the fractional quasiparticles in the 
bulk satisfy the so-called “semionic” statistics, as expected 
for Laughlin u = 1/2 state. In the topological quantum held 
theory, quasiparticle statistics is encoded in the modular ma¬ 
trices which describe the action of modular transformation 
on the topological groundstates"^^"^^. The 5—matrix contains 
the mutual statistics information of the anyonic quasiparticles, 
such as quantum dimensions and fusion rules between differ- 


III. DISCUSSION AND CONCLUSION 

We have established a phase diagram for the Bose-Hubbard 
model for kagome system at 1/3 Ailing number with the emer¬ 
gent FQH phase in the intermediate regime. We also address 
the nature of the quantum phase transitions between the FQH 
phase and solid phases (see Appendix D). We utilize several 
quantities, such as entanglement entropy, correlation length 
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and groundstate wavefunction fidelity. The numerical evi¬ 
dences signal the first order character of the phase transition 
between FQH phase and the strip phase as well as the charge 
density wave state. Moreover, we also find that in all three 
phases, the obtained correlation length is much smaller than 
the cylinder width Ly, which confirms that our DMRG calcu¬ 
lation offers a reliable phase diagram for the thermodynamic 
limit. 

Regarding the laboratory realization for the emergent FQH 
state, a natural experimental setting for Bose-Hubbard physics 
is ultracold atomic gases^^. One advantange of our model 
is that it only contains real nearest-neighbor hoppings and 
density-density interactions. Moreover, we also note that 
there are several existing candidates of spin—1/2 materi¬ 
als with kagome structure, such as BaCusV 20 s{ 0 H) 2 ^^, 
Cus{Mg, Zn)07{0H)2 • H20^\ Cu^V207{0H)2 • 
Rh 2 Cu 2 ,SnFi 2 ^^ and DysRu 4 Ali 2 ^^, and each of them has 
its own interactions deserving to be studies more carefully un¬ 
der magnetic field for possible detecting of the exotic mag¬ 
netization plateaus. But the fine tuning of second and third 
nearest neighbor interactions in such materials may be diffi¬ 
cult to achieve. So the material-realization of our proposed 
FQH phase in condensed matter setting may depend on syn¬ 
thesising more kagome materials in the future. 

In conclusion, we have presented a global phase diagram 
of an extended Bose-Hubbard model on the kagome lattice at 
fractional one-third filling. Importantly, the interplay between 
the underlying lattice and strong interaction gives birth to a 
fractional quantum Hall (FQH) liquid phase, even though the 
non-interacting band structure is topological trivial. The FQH 
phase we present here provides a “proof of the principle” ex¬ 
ample of interaction-driven quantum anomalous Hall effect 
with time-reversal symmetry (TRS) spontaneously breaking. 
We also provide complete characterization of the universal 
properties of the FQH phase, including ground state degen¬ 
eracy, topological entanglement spectrum, fractionally quan¬ 
tized Chern number and anyonic quasiparticle statistics. To 
our best knowledge, this is the first example of a TRS break¬ 
ing topological phase at one-third filling on kagome lattice 
system. 

We believe our current work will inspire upcoming research 
efforts both in theoritical and experimental fields. From the 
theoretical side, the present calculations show that the model 
of strongly interacting hard-core bosons can harbor rich and 
interesting phases through interaction engineering. It is also 
interesting to study the interacting particles to be spinless or 
spinful fermions, since the possible topological liquid phase in 
fermionic models has been sought for a long time^^’^^ From 
the experimental side, our work will provide direction and in¬ 
sight in searching for the topological liquid phase in realis¬ 
tic materials with kagome lattice structure, or by engineering 
such systems in ultracold atomic settings. Thus our present 
findings would provide both theoristes and experimentalists 
a rich playground in searching of new topological phases in¬ 
duced by strong interaction. 
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Appendix A: Method 

In this paper, the calculations are based on the density- 
matrix renomralization group (DMRG) algorithm on cylinder 
geometry^^’^"^ and the exact diagonalization (ED) on torus ge¬ 
ometry, both of which have been proven to be effective and 
complementary tools for studying realistic models containing 
arbitrary strong and frustrated interactions. On one hand, ED 
is staightforward in identifying the groundstate degeneracy on 
compactified spaces. But the drawback is that with the expo¬ 
nential growing of the Hilbert space, the accesible systems 
are limited to smaller sizes, up to Ng = 36 for this study. 
On the other hand, DMRG calcualtion allows us to obtain ac¬ 
curate groundstates and related entanglement measurements 
on much larger system sizes beyond the ED limit. More¬ 
over, DMRG calculation also has the advantages of probing 
ground states with spontaneous symmetry breaking and topo¬ 
logical ordering. The DMRG calculations on long cylinder 
tend to automatically select the groundstates with minimal 
entropy^^, which is helpful to study the fractionalized quasi¬ 
particle statistics in topological ordered states^ 


1. Details of DMRG Calculation 

We study the cylinder system with open boundaries in the 
X direction and periodic boundary condition in the y direc¬ 
tion. The available system sizes are cylinders of circumfer¬ 
ence Ly = 3,4, 5,6 (in unit of unit cell). Eor the largest sys¬ 
tem width {Ly = 6), we keep up to M = 8400 U{1) states 
and reach the DMRG truncation error 5 x 10“^. 

The entanglement entropy and spectrum can be easily ob¬ 
tained in the DMRG. By partitioning the system into sub¬ 
systems A and B, the groundstate wavefunction \ijj) can 
be decomposed according to Schmidt decomposition j?/) = 

-I /Q 

where are eigenvalues of the reduced 
density matrix pA of subsystem A. Thus the entanglement en¬ 
tropy can be defined as Sa = —In Pa] = — 

The eigenvalues log{\i} plotted against the relative momen¬ 
tum quantum number /S.ky of the subsystem A, is defined as 
the entanglement spectrum^^. 


2. Adiabatic DMRG and Fractionally Quantized Chern 
Number 

We have used the numerical fiux insertion experiment based 
on the adiabatical DMRG simulation to detect the topological 
Chern number of the bulk system^^’^^’"^^. To simulate the fiux 0 
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threading in the hole of a cylinder, we impose the twist bound¬ 
ary conditions along the y direction with replacing terms 
bl,br + h.c. e'^^^'^bl,br + h.c. for all neighboring (r, r') 
bonds with hoppings crossing the y-boundary in the Hamilto¬ 
nian (Eq. 1). The charge pumping from one edge to the other 
edge can be computed from {AQ{0)) = Tr[pL{0)Q{0)], 
where Q{0) is the U{1) quantum number and Pl{0) is re¬ 
duced density matrix of left half system. Due to the quantized 
Hall response, the Chern number of ground state is equal to 
the charge pumping by threading a 6> = 27r flux^^. To realize 
the adiabatic flux insertion, we use the step of flux insertion 
as AO = 0.257r. 


Appendix B: Classical Phase Diagram 

In this section, we discuss the phase diagram of the model 
without hopping term t = 0, where the system reduces 
to a classical Ising model with competing antiferromagnetic 
interactions^^. We compare the energy of four states, i.e., the 
q = (0,0) state, the q = (0, tt) state, the stripe state, and the 
\/3 X ^/3 state. The pattern of these states are shown in the 
inset of Fig. 6(a). Interestingly, we And that when both V 2 , V 3 
are larger than Vi, the stripe state always has the lowest en¬ 
ergy; otherwise, the q = (0, 0) and q = (0, tt) states compete 
depending on the strengths of V 2 and V 3 (see Fig. 6(a)). 
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FIG. 6: Classical phase diagram in {14,1^2,14} parameter space. 

(a) The phase diagram plotted in V2/V1 and V3/V1 parameter space. 
The insets show the unit cells of the different phases, where the 
large and small circles denote the occupied and the unoccupied sites. 

(b) The phase diagram plotted in Vi and V 2 = Vs parameter space. 

It is interesting to compare this classical phase diagram with 
the quantum phase diagram as shown in Fig. 1(a). To compare 
them, we set V 2 = V 3 as shown in Fig. 6 (b), where the system 
shows two phase regions, a stripe phase and a charge den¬ 
sity wave phase region with the boundary at ^2 = ^3 = ^!. 
Interestingly, the charge density wave states with q = (0,0) 
and q = (0,7r) are degenerated in the case of < V 2 - In the 
quantum case, the degeneracy is lifted with the q = ( 0 , 0 ) state 
having the lower energy. The FQH phase obtained by DMRG 
appears in the transition region between the stripe phase and 
the q = (0,0) phase of the classical phase diagram. The 
emergence of the topological phase appears seems to arise as 
a result that quantum fluctuations destory long-ranged orders 
around the transition region. This could serve as a guiding 
principle for flnding topological phases in other models and/or 
on other lattices. 



FIG. 7: Entanglement entropy scaling with Ly. 


Appendix C: Topological entanglement entropy 

For a gapped quantum phase with topological order, 
the topological entanglement entropy (TEE) 7 is proposed 
to characterize the non-local entanglement^^’^^. Generally 
speaking, entropy has the form S = aLy — 7 , where Ly is the 
boundary of the subsystem, and a is a non-universal constant. 
While a positive 7 is a correction to the area law of entangle¬ 
ment and reaches a universal value determined by total quan¬ 
tum dimension V of quasiparticle excitations as 7 = In D. For 
the z/ = 1/2 Faughlin state, the quantum dimension of each 
quasiparticle is 1 (see main text), leading to the total quantum 
dimension V = \/2 and thus the TEE 7 = In D = ^ In 2. 

By using the DMRG simulations, we obtain the minimal 
entropy state^^’^^ with spontaneously broken time-reversal 
symmetry and calculate the corresponding Von Neuman en¬ 
tanglement entropy. The converged entropy are available for 
Ly = 3,4, 5,6 cylinders. For Vi = 0.5, V2 = 1-0, V3 = 1.0, 
we make a linear fitting of the entropy data for = 3 ^ 4 , 5,6 
cylinders, and And the TEE 7 « 0.343 ± 0.075. (If we make 
the linear fitting based on data for Ly = 4, 5, 6 (not shown), 
the obtained result is 7 ^ 0.241 ± 0.098). Despite some un¬ 
certainty in the fitting, the obtained TEE approaches the pre¬ 
diction of the 1 / = 1/2 Laughlin state 7 = ^ In 2 « 0.346. 

Appendix D: Quantum phase transitions 

In order to uncover the nature of correponding phase transi¬ 
tions between FQH phase and solid phases, we inspect several 
quantities that are expected to be sensitive to a phase tran¬ 
sition, such as the entanglement entropy S, and correlation 
length Both quantities are expected to show a finite jump 
when crossing a first order transition. We also calculate the 
groundstate wavefunction fidelity F = \{'ip{V)\'ip{V 
(V is some parameter in Hamiltonian)^"^, which can faithfully 
describe the first-order transition or energy level crossing. 

We show the results along the reference line in Fig. 1 (in 
the main text), by fixing Vi = 0.5 and varying V 2 = Vs- In 
Fig. 8 (a), it is found that the entanglement entropy shows 
a sharp jump around V 2 = V 3 ^ 0.8 and a drop around 
V 2 = V 3 ~ 1.5. We also observe the similar behavior when 
looking at the correlation length, as shown in Fig. 8 (b). Both 
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FIG. 8: (a) Entanglement entropy and (b) correlation length versus 
1 / 2 (= V 3 ) by setting Vi = 0.5. The calculations are performed on 
Ly — A cylinder, by keeping M = 800 (brown dots) and 1600 (blue 
dots) states, (c) Wavefunction fidelity plotted as y 2 (= Vs) by setting 
Vi = 0.5. The calculations are performed on = 4 cylinder, by 
keeping M = 800. 


of these two measurements signals a direct first order phase 
transition between FQH phase and solid phases. In addition, 
we also find that, between the FQH phase and charge den¬ 
sity wave q = (0,0) phase, there exists a narrow window 
for the coexistence of both FQH nature and the charged or¬ 
der, as shown by the shaded area in the phase diagram. We 
have checked that, in this intermediate regime, the ground- 
state hosts the quantized Chern number C = 1/2, but devel¬ 
ops weak charge order. From the wavefunction fidelity in Fig. 
8(c), it is shown a first-order transition between the coexis¬ 
tence region and the FQH phase (q = (0, 0) phase). 


Moreover, when studying the topological order on cylin¬ 
der geometry with finite width Ly, the correlation length ^ of 
the ground state offers a natural consistency check for the as¬ 
sumption that the value of Ly is large enough to be represen¬ 
tative of the thermodynamic limit. The correlation length is 
defined by ^ = — In |ei/e 2 |, where ei ^2 are two largest eigen¬ 
value from transfer matrix^^. If Ly is much larger than we 
indeed expect finite size effects to be very small. Indeed, in 
Fig. 8(b), the condition ^ is satisfied so that our DMRG 
calculation offers a reliable and relevant results for thermo¬ 
dynamic limit. Based on these measurements, we expect the 
finite size effect should be small in our calculations. 
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